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Summary. In a recent paper [17] we conjectured that for ferromagnetic Heisen¬ 
berg models the smallest eigenvalues in the invariant subspaces of fixed total spin 
are monotone decreasing as a function of the total spin and called this property fer¬ 
romagnetic ordering of energy levels (FOEL). We have proved this conjecture for the 
Heisenberg model with arbitrary spins and coupling constants on a chain [17, 20]. 
In this paper we give a pedagogical introduction to this result and also discuss some 
extensions and implications. The latter include the property that the relaxation time 
of symmetric simple exclusion processes on a graph for which FOEL can be proved, 
equals the relaxation time of a random walk on the same graph. This equality of 
relaxation times is known as Aldous’ Conjecture. 


1 Introduction 

The ferromagnetic Heisenberg model is the primordial quantum spin model. It has 
been studied almost continuously since it was introduced by Heisenberg inl926. In 
the course of its long history, this model has inspired an amazing variety of new 
developments in both mathematics and physics. The Heisenberg Hamiltonian is 
one of the basic, non-trivial quantum many-body operators, and understanding its 
spectrum has been a guiding problem of mathematical physics for generations. 

A lot of attention has been given to the Bethe-Ansatz solvable one-dimensional 
spin-1/2 model, which has an infinite-dimensional algebra of symmetries [9]. The 
results we will discuss here are not related to exact solutions but there is an essential 
connection with the 51/(2) symmetry of the model, much in the spirit of the famous 
result by Lieb and Mattis ([15], see also [14, footnote 6]). The Lieb-Mattis Theorem 
proves “ordering of energy levels” for a large class of antiferromagnetic Heisenberg 
models on bipartite lattices. Namely, if the two sublattices are A and B, and all 
interactions within A and B are ferromagnetic while interactions in between A 
and B are antiferromagnetic, then the unique ground state multiplet has total spin 
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equal to |<Sa — Sb\, where Sa and Sb are the maximum total spins on the two 
sublattices. Moreover, the minimum energy in the invariant subspace of total spin 
S, for S > |iSa — <Ss|, is monotone increasing as a function of S. The most important 
example where this theorem provides useful information is the usual antiferromagnet 
on a bipartite lattice with equal-size sublattices. Then the ground state is a unique 
spin singlet, and the minimum energy levels for each possible total spin S, are 
monotone increasing in S. Our aim is a similar result for ferromagnets. To be able 
to state the ferromagnetic ordering of energy levels (FOEL) property precisely, we 
hrst give some dehnitions. 

Let /I be a finite connected graph with a set of vertices or sites, x, that we will 
also denote by A and a set E of unoriented edges, or bonds, {xy). We will often write 
X ^ y £ A to signify that the edge (xy) is present in A. In many physical examples 
one has A C 

Each site x £ A has a quantum spin of magnitude Sx £ {1/2,1, 3/2,...}, asso¬ 
ciated with it. The state space at x is 2sx + 1-dimensional and we denote by S'{, 
i = 1,2,3, the standard spin-s^, matrices acting on the a:th tensor factor in the 
Hilbert space Ti. — ®The isotropic (also called XXX) ferromagnetic 
Heisenberg Hamiltonian on A is given by 

Ha = - Y. J-yS.-Sy, ( 1 ) 

where the real numbers Jxy are the coupling constants, which we will always as¬ 
sume to be strictly positive (that they are positive is what it means to have the 
ferromagnetic Heisenberg model). This model is widely used to describe ferromag¬ 
netism at the microscopic level whenever itinerant electron effects can be ignored. 
Examples are magnetic domain walls and their properties and a variety of dynamical 
phenomena. 

The spin matrices generate an irreducible representation of SU (2) at each vertex. 
This representation is conventionally denoted by An important feature of the 

Hamiltonian (1) is that it commutes with SU{2) via the representation 

(g) (2) 

x€A 

or, equivalently, with the total spin matrices defined by 

S*a = Y^-’ i = 1,2,3. 

x^A 

and hence also with the Casimir operator given by 


C = Sa-Sa. 

The eigenvalues of C are S{S -I- 1), S' = Smin, Smin -I- 1,..., Smax = Yjx£A which 
are the spin labels of the irreducible representations that occur in the direct sum 
decomposition of the tensor product representation (2) into irreducible components. 
The value of Smin is usually 0 or 1/2, but may be larger if one of the Sx is greater 
than Smax/2. The decomposition into irreducible components can be obtained by 
repeated application of the Clebsch-Gordan series; 

^(»l) (g| ^ ®2|) 0 £)(|si —*2 l-l-l) ... 0 £)(®l+®2) 


( 3 ) 
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The label S is called the total spin, and the eigenvectors of the eigenvalue S{S + 1) 
of C, are said to have total spin S. Let denote the corresponding eigenspace. 
Since C commutes with Ha, the spaces are invariant subspaces for Ha- For 
any hermitian matrix H leaving the spaces invariant we define 

E{H,S) = min spec 

By Ferromagnetic Ordering of Energy Levels (FOEL) we mean the property 

E{H,S) < E{H,S'), if S'' < S. 

for all S and S' in the range [Smin, Smax]- 

In particular, if Ha has the FOEL property it follows that its ground state 
energy is E{Ha, Smax), which is indeed well-known to be the case for the Heisenberg 
ferromagnets. Moreover, since the multiplet of maximal spin is unique, FOEL also 
implies that the gap above the ground state is E{Ha, Smax — 1) — E{Ha, Smax), 
which is well-known for translation invariant Heisenberg ferromagnets on Euclidean 
lattices. 

Conjecture 1. All ferromagnetic Heisenberg models have the FOEL property. 

The FOEL property and the Lieb-Mattis theorem applied to a spin-1 chain of 5 
sites is illustrated in Figure 1. 

Our main result is a proof of this conjecture for the special case of arbitrary 
ferromagnetic Heisenberg models on chains, i.e., one-dimensional model [17, 20, 18]. 

Theorem 1. FOEL holds for ferromagnetic XXX spin chains, i.e., for all 

^--1 1 

H = -y J.,.+i(-5,. • S.+i - 1), (4) 

—' Sa;Sa, + l 

X = 1 

for any choice of Sj, G {1/2, 1, 3/2,...} and Jx,x +1 > 0. 


2 Proof of the Main Result 

Our proof of (Theorem 1) proceeds by a finite induction argument for a sequence 
of models with Hamiltonians Hk ~ Hf., 1 < k < N, on Hilbert spaces Hk, with the 
following properties: 

(i) There is a unitary representation of SU{2), Uk, on Hk, that commutes with 
Hk. 

(ii) There are isometries 14 : Hk+i —> Hk ®C^, interwining the representations 

Uk+i and Uk ® i.e., VkUk+i{g) = {Uk{g) ® D^^^^\g))Vk, for all g G SU{2), 

and such that 

Hk+i > V:{Hk ® 11)14 

(iii) Hi has the FOEL property. 

(iv) For every S' = 0,1/2, 1, 3/2, ..., for which 77^®^ / {0}, we have 


E{Hk+i,S+1/2) < E{Hk,S). 
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Spectrum of SU(2) - symmetric Heisenberg model : j = 1, L = 5 





- - - - - - - - --• 4 
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S*^*-eigenvalue / M-sector 

Fig. 1. The spectrum of a ferromagnetic Heisenberg chain consisting of 5 spin-1 
spins, and with constant couplings. On the horizontal axis we have plottted the 
eigenvalue of the third component of the total spin. The spectrum is off-set so that 
the ground state energy vanishes. The arrows on the right, with label S, indicate 
the multiplets of eigenvalues E{H,S), i.e., the smallest eigenvalue in the subspace 
of total spin S. The monotone ordering of the spin labels is the FOEL property. On 
the left, we have indicated the largest eigenvalues for each value of the total spin. 
The monotone ordering of their labels in the range 1,..., 5, is the content of the 
Lieb-Mattis theorem applied to this system. 


We will first present the induction argument using the assumptions (i)-(iv), and 
then construct the sequence Hk satisfying these four assumptions. This argument is a 
generalization of results in [12] and [17]. The sequence of Hamiltonians will, roughly 
speaking, be a sequence of systems of increasing size, starting with the trivial system 
of a single spin. Property (i) simply means that all models will have isotropic inter¬ 
actions. Property (ii) will closely guide the construction of our sequence. Property 
(iii) will be trivial in practice, since Hi will be a multiple of the identity on Tli in 
our applications. Property (iv) has a nice physical interpretation at least in some 
of the examples we will consider (see Section 5.2). It is our (in)ability to prove (iv) 
that limits the range of models for which we can prove FOEL. 

Theorem 2. Let {Hk)i<k<N, be a sequence of Hamiltonians satisfying properties 
(i)-(iv). Then, for all k, Hk has the FOEL property. 
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Proof. Since Hi has the FOEL property by assumption, it is sufficient to prove the 
induction step. Consider the following diagram: 

E{Hk,S) >1 E{Hk,S + l) > E{Hk,S + 2) 

E{Hk+i,S + ^) >4 E{Hk+i,S+^) > .B(Hfc+i,S + |) 

The inequality labeled 1 is FOEL for Hk, and inequality 2 is property (iv) assumed 
in the theorem. We will prove inequality 3 (using inequality 1) and, combined with 
inequality 2 this implies inequality 4, which is the induction step. 

As before, we use superscripts to Hilbert spaces to denote their subspaces of 
fixed total spin. To prove inequality 3, we start from the variational principle: 


E{Hk+i,S+ 1/2) = inf (<(>,Lffc+i<(>) 


> inf {/},V:{Hk(St 2 )Vk/}) 


> 


inf 


{i>, (Hk ® ]l2)V>) 


The first inequality uses the fact that 14 is an isometry and property (ii). For the 
second inequality we enlarged the subspace over which the infimum is taken. 

Now, we use the Clebsch-Gordan series (3) to see that {Hk <8) C 

(Ltf ’ © Ltf© C^. Therefore 

E{Hk+i, S + 1/2) > mm{E{Hk, S), E{Hk, S + l)} = E{Hk, S + 1). 

Clearly, Hk © 12 restricted to 0 has the same spectrum as Hk 

(S) rs+1') 

restricted to H], ' © H], '. The last equality then follows from inequality 1, i.e., 

the induction hypothesis. This concludes the proof of Theorem 2 


For the proof of Theorem 1 we will apply Theorem 2 to the sequence {Hk)i<k<N, 
with N = 2S'max, Hi — 0, and Hn = H, constructed as follows: for each k = 
1,..., Af — 1, the model with Hamiltonian Hk+i is obtained from Hk in one of two 
ways: either a new spin 1/2 is added to right of the chain, or the magnitude of 
the rightmost spin is increased by 1/2. In both cases, Smax goes up by 1/2 at each 
step, hence N = s^. Each Hk is of the form (4), and we have written the 

interactions in such a way that the coupling constants Jx,x+i can be taken to be 
independent of k, although this is not crucial since all arguments work for any choice 
of positive coupling constants at each step. The parameters that change with k are 
thus L and the set of spin magnitudes {sx)x=i. To be explicit, the two possible ways 
of deriving Hk+i from Hk are summarized in Table 2. 

The Hamiltonians are of the form 


Hk — ^ ) Jx,x.\-1 

X = 1 


1 


Sx{k)sx+i{k) 


Sx ' *S'x + l 1 


(5) 


where S'),, i = 1, 2, 3, are the 2sx{k) + 1 dimensional spin matrices. To simplify the 
notation, the dependence on k will often be omitted further on. 

We now have a uniquely defined sequence of Hamiltonians {Hk)i<i<N, with 
Hi = 0 and Hn = H. Next, we proceed to proving the properties (i)-(iv). Property 
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(i) is obvious by construction. Property (iii) is trivial since Hi = 0. To verify property 

(ii) , we need to distinguish the two cases for the relation between Hk and Hk+i, as 
given in Table 2. 

For Case I, Uk+i = and we can take the identity map for V. Property 

(ii) follows from the positivity of the additional interaction term in Hk+i'- 

Hk+i = Hk + Jli,,Li,+i (- ■ Slu+i — 1 

For Case II, we have Hk = Hi <3 and Hk+i = Hi ® for 

some I < k, possibly I = 0,Ho = C. Since = Sl^ + 1/2, there is a (up to a 

phase) unique SU{2) intertwining isometry W : ^ C)C^, namely 

the W that identifies the spin subrepresentation in From the 

intertwining property, the irreducibility of the spin representations, and the SU{2) 
commutation relations one deduces that there is a constant c such that 


W*(Si^ (fc) ® fl)VF = cSi, (k + l), i = 1, 2, 3. 


The constant c is most easily determined by calculating the left and right hand sides 
on a highest weight vector (a simultaneous eigenvector of C and with eigenvalues 
S{S + 1) and S, respectively). One hnds 

SL^jk) 

SLk{k + 1 ) 


Now, take V = Ih, ® W. It is then straightforward to check that 


V* 


sli._iSli. {k) 


Sl 


■Sl 


112 1 ^ = 


1 


SLt,-lSL^{k + 1 ) 


Sl,_, ■ Sr 


where the spin matrices on the left hand side are of the magnitude determined by 
SLk-i and sl^ (k) , while on the right hand side they are the magnitudes of the spins 
are sl^-i and sz,;,(fc + l). 

To prove Property (iv), we start by observing that 


spec{Hk\^is)) = spec(//fc|^(s)) 

where V/ ' is the subspace of Hk of all highest weight vectors of weight S. This is 
an invariant subspace for Hk and for every eigenvalue of Hk\^(s) there is at least 

one eigenvector in Let d{k,S) denote the dimension of 

Property (iv) will be obtained as a consequence of the following proposition and 
a version of the Perron-Frobenius Theorem. 


Table 1. Summary of the fe-dependence of the sequence of models used in the proof 
by induction of Theorem 2. 


parameter 

Case I 

Case II 

L 

Lk+l ~ Lk + i 

Lfc+i = Lk 


SLk+i{k) = 0, SLk+i{k + 1) = 1/2 

SLk (fc + 1) = SLk (k) + 1/2 

H 

Hk+i ^Hk(E)C^ 

Lffc+i = I/(Lffc ® C") 
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Proposition 1. We have d{k + 1,5'+ 1/2) > d{k,S) and there are bases for 
such that the matrices of Hk\^{s) with respect to these bases have the 

following properties: 

^(k.s) ^ forl<i^j < d{k, S),l<k<N 
^(k+i,s+i/ 2 ) ^ l<i,j < d{k, S),l<k<N-l. 

For reasons of pedagogy and length, we will give the complete proof of this 
proposition only for the spin 1/2 chain. The proposition provides the assumptions 
needed to apply a slightly extended Perron-Frobenius theorem (see, e.g., [23]), which 
we state below. 

The standard Perron-Frobenius Theorem makes several statements about square 
matrices with all entries non-negative, which we will call a non-negative matrix 
for short. Recall that a non-negative matrix A is called irreducible if there exists 
an integer n > 1 such that the matrix elements of A" are all strictly positive. 
The standard results are the following: (i) every non-negative matrix has a non¬ 
negative eigenvalue equal to its spectral radius (hence it has maximal absolute value 
among all eigenvalues), and there is a corresponding non-negative eigenvector (i.e., 
with all components non-negative); (ii) if A is an irreducible non-negative matrix 
there is a unique eigenvalue with absolute value equal to the spectral radius of A, 
which is strictly positive and has algebraic (and hence geometric) multiplicity 1. Its 
corresponding eigenvector can be chosen to have all strictly positive components. 

If A is a square matrix A with all off-diagonal matrix elements non-positive, we 
will call A irreducible if there exists a constant c such that cl — A is irreducible 
according to the previous definition. From the standard Perron-Frobenius Theorem 
it immediately follows that the eigenvalue with smallest real part of an irreducible 
matrix in the last sense is real, has algebraic (and hence geometric) multiplicity 1, 
and that the corresponding eigenvector can be chosen to have all components strictly 
positive. In the following, we will repeatedly use the information provided by the 
standard Perron-Frobenius Theorem as described above without further reference. 
Let specrad(A) denote the spectral radius of a square matrix A. 

Lemma 1. Let A = {oij) and B = {bij) be non-negative nxn matrices, and assume 
that Oij < bij, for all 1 < i, j < n. Then 

specrad(A) < specrad(i3). (6) 

If B is irreducible and there is at least one pair ij such that Oij < bij, then 

specrad(A) < specrad(i3). (7) 

Since the spectral radii are also the eigenvalues of maximal absolute value, the same 
relations holds for these eigenvalues. 

Proof. Let r = specrad(A). Then A has a non-negative eigenvector, say v, with 
eigenvalue r. If Oij < bij, for all 1 < *,/ < n, it is clear that there is a non-negative 
vector w such that 

Bv = rv w. (8) 

This relation implies that ||H^|| > r*', for all positive integers k and, hence, 
specrad(i3) > r. This proves (6). 
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To prove (7) for irreducible B such that atj < bij for at least one pair of indices, 
let fc be a positive integer such that is strictly positive. This implies that B'^v 
has all strictly positive components. From this it is easy to see that the non-negative 
w such that 

S k k I 

V = r V + w 

cannot be the zero vector. Therefore there is 2 £ R with all strictly positive compo¬ 
nents such that 

Dfc +1 fc +1 . 

B v = r v + z. 

Since 2 is strictly positive, there exists e > 0 such that, ev < z componentwise, and 
therefore we can find <5 > 0 such that 

B'‘+\ = (r + + z', 

with 2 ' non-negative. We conclude that specrad(B) > r + S > specrad( 2 l). 

Note that the argument that proves this lemma could also be used to give a 
lower bound for the difference of the spectral radii. Since we do not need it, we will 
not pursue this here. The next theorem is an extension of Lemma 1. 


Theorem 3. Let A = (aij) and B = {bij) be two square matrices of size n and m, 
respectively, with n < m, both with all off-diagonal matrix elements non-positive, 
and sueh that bij < Oij, for 1 < *, j < n. Then 

inf spec(B) < inf spec(>l) (9) 

If B is irredueible and either (i) there exists at least one pair ij, 1 < i,j < n, sueh 
that bij < aij; or (ii) bij < 0, for at least one pair ij with at least one of the indices 
i or j > n, then 

infspec(B) < infspec(yl). (10) 


Proof. Let c > 0 be a constant such that the matrices A' = (aT) = clln — A and 
B' = {b'ij) = dm — B are non-negative. Define A" to be the mx m matrix obtained 
by extending A' with zeros: 


(a") = A" 


'A' o' 
0 0 


It is easy to see that a'ij < b'ij, for 1 < i, j < m. Therefore, we can apply Lemma 1 
with A" playing the role of A, and B' playing the role of B. Clearly, specrad(T") = 
c — inf spec( 2 l) and specrad(B') = c — inf spec(B). Therefore, this proves (9). 
Similarly, (10) follows from the additional assumptions and (7). 


Proof of Theorem 1: The remaining point was to prove property (iv) needed in 
the assumptions of Theorem 2. We use Proposition 1, which we will prove in the 
next section, and apply Theorem 3 with A = and B = This 

completes the proof. 



Ordering of Energy Levels in Heisenberg Models and Applications 9 

3 The Temperley-Lieb basis. Proof of Proposition 1. 

In the proof of Theorem 1 in the previous section we used the matrix representation 
of the Hamiltonians restricted to the highest weight spaces given by Proposition 
1. We now give the complete proof of that proposition for the spin 1/2 chain and 
sketch the proof in the general case. 

The main issue is to hnd a basis of the highest weight spaces with the desired 
properties. Fortunately for us, such a basis has already been constructed and we 
only need to show that it indeed had the properties claimed in Proposition 1. For 
the spin 1/2 chain we will use the Temperley-Lieb basis [22], and for the general 
case its generalization to arbitrary spin representations introduced by Frenkel and 
Khovanov [6]. 

3.1 The basis for spin 1/2 

We start with the spin 1/2 chain, i.e., Sx = 1/2, for all x. In this case Smax = k/2 
and is the subspace of consisting of all vectors ip such that S^ip = Sip 

and S'^ip = 0. Let n be the “spin-deviation” dehned as S = k/2 — n. Then, n is 
a non-negative integer. The case n = 0 is trivial since dimV^*^^^ = 1, namely just 
the mutiples of the vector |-|-) ® |-|-) (g) • • • (g) |-|-), where |±) is the basis of that 
diagonalizes S^. For n > 1, the basis vectors are a tensor product of n singlet vectors 
^ = |-l-) <8 | —) — |—) ® |-l-), accounting for two sites each, and k — 2n factors equal to 
|-|-). Such vectors are sometimes called Hulthen brackets. It is clear that any such 
factor is a highest weight vector of weight k/2 — n, just calculate the action of S'® 
and S®" on such a vector. They are not linearly independent however, except in 
the trivial case k = 2. The contribution of Temperley and Lieb was to show how 
to select a complete and linearly independent subset, i.e., a basis. How to select 
the Temperley-Lieb basis, is most easily explained by representing the vectors by 
configurations of n arcs on the k vertices 1,..., n. The arcs are drawn above the 
line of vertices as shown in Figure 2. Each arc represents a spin singlet and each 
unpaired vertex represents a factor |-|-). The vectors (conhgurations of arcs) selected 
for the basis are those that satisfy two properties: (i) the arcs are non-crossing, (ii) 
no arc spans an unpaired vertex. The resulting set is a (non-orthogonal) basis. E.g., 
the basis for k = 5 and n = 2 is shown in Figure 2. We will use, a,l3,.. ., to denote 
arc configurations that obey these rules, and by the corresponding basis vectors will 
be denoted by ja) ,[/?),.... We will use the notation \xy] € a to denote that the arc 
connecting x and y is present in a. 

Proof of Proposition 1 for the spin 1/2 chain. The action of the Hamiltonian 
on the basis vectors has an appealing graphical representation. We can write the 
Hamiltonian as 

fc-i 

Pkk — 2 ^ ] Tx.x + l + l 

X = \ 

where Ux,x+i = —C ® C* which, up to a factor —2, is the orthogonal projection onto 
the singlet vector acting on the xth and x + 1st factor in the tensor product. The 
Ux,x+i form a representation of the Temperley-Lieb algebra with parameter q = 1 
(see, e.g, [10]). It is a straightforward calculation to verify the action of Ux,x+i on a 
basis vector ]a): (i) if both x and a; -I- 1 are unpaired vertices in a, Ux,x+i ja) = 0; 
(ii) if [x,x -I- 1] G a, we have Ux,x+i |«) = —2 ]a); (hi) if [uw] G a, with exactly one 
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Fig. 2. The possible configurations of 2 arcs on 5 vertices. 


of the vertices u and v equal to x or a; + 1, we have Ux,x+i la) = \P), where /? is 
obtained form a by removing [uu] and adding [x,x + 1]; (iv) if \ux\ and [a; + l,w] 
are both present in a, we have Ux,x+i |a) = |/?), where (5 is obtained form a by 
removing [ua:] and x + l,u], and adding \uv\ and [a;,a: + 1]. 

The action of Ux,x+i on the vector |a) can be graphically represented by placing 
the diagram shown in Figure 3 under the diagram for a, and read off the result using 
the graphical representation of the rules (i)-(iv) shown in Figure 4. The action of 
the Hamiltonian is then obtained by summing over x as shown in Figures 5 and 6 . 



Fig. 3. The graphical representation of Ux,x+i. 


The important observation is the action of the Hamiltonian on a basis vector 
I a) yields a linear combination of basis vectors with non-positive coefficients except 
possibly for the coefficient of [a) itself, which has the opposite sign resulting from 
the “bubble” in the graphical representation. This means that all off-diagonal matrix 
elements are non-positive as claimed for the matrices in the proposition. 

The second will follow from the observation that is a submatrix of 

j^k+i,s+i/ 2 ^ Note that the spin deviation for and is the same, say 

n. Let us order the basis elements of so that all a where the last vertex, 

fc -I- 1, is unpaired, are listed first, and consider the afi matrix element of Hk+i 
for such a and /?. Then, it is easy to see that there are no contributions from the 
k,k+l term in Hamiltonian, since its action results in non-zero coefficients only for 
configurations where fc -|-1 belongs to an arc. This means that these matrix elements 
are identical to those computed for Hk for basis vectors labeled a' and fi' obtained 
from a and (3 by dropping the last vertex, fc -|- 1 which is unpaired. 
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Fig. 4. The graphical rules (i)-(iv) for 
basis vector. 



Fig. 5. Action of the Hamiltonian of the 
alized Hulthen bracket, for L = A, k = 1. 


the action of Uj;^x+i on a Temperley-Lieb 



spin-1/2 XXX or XXZ chain on a gener- 


This completes the proof of Proposition 1 in the case of the pure spin 1/2 chain. 

Q.E.D. 

3.2 The basis for higher spin 

We are looking for a basis of the space of highest weight vectors of weight S of 
the spin chain with Hilbert space Hk- Equivalently, we may look for a basis of the 
SU(2) intertwiners ^ TLk- There is a graphical algebra of such intertwiners 
with a very convenient basis, the dual canonical basis, introduced by Frenkel and 
Khovanov [6]. This is the basis we will use, but we will present it as a basis for the 
subspaces V/ ’ of highest weight vectors. 
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Fig. 6. Action of the Hamiltonian of the spin-1/2 XXX or XXZ chain on a gener¬ 
alized Hulthen bracket, for k = 6, n = 2. 


The state space at site x can be thought of as the symmetric part of 2sx spins-i. 
We can label the 2sx + 1 states by the Ising conhgurations 

ITT ••• T>, lit ••• T>, lUT ■■■ T), •••, • 

where each configuration stands for the equivalence class up to re-ordering of all 
configurations with the same number of down spins. E.g., [TTTTT) is the vector 
normally labelled as m) = |5/2,1/2), and not the tensor |J,) ® |J,) (g) ||) ® |t) ® |t}. 
The states for a chain of L spins of magnitudes si,..., si, are then tensor products 
of these configurations. We shall call such vectors ordered Ising eonfigurations. These 
tensor product vectors, in general, are not eigenvectors of the Casimir operator S, 
i.e., they are not of definite total spin. Suitable linear combinations that do have 
definite total spin are obtained by extending the Hulthen bracket idea to arbitrary 
spin as follows. Start from any ordered Ising configuration such that 2M = # T 
i- Then, look for the leftmost J, that has a T to its left, and draw an arc 
connecting this j to the rightmost |, left of it. At this point, one may ignore the 
paired spins, and repeat the procedure until there is no remaining unpaired j, with 
an unpaired / to its left. This procedure guarantees that no arcs will cross and 
no arc will span an unpaired spin. The result, when ignoring all paired spins, is 
an ordered Ising conhguration of a single spin. See Figure 7 for an example of this 
procedure. The result is a basis for the spin chain consisting entirely of simultaneous 
eigenvectors of the total spin and its third component, with eigenvalues S and M, 
respectively. The value of M is 1/2 times the difference between the number of up 
spins and the number of down spins in the ordered Ising configuration. The total- 
spin S is equal to S minus the number of pairs. Clearly, the highest weight vectors 
are then those that have no unpaired j, i.e., the ordered Ising configuration consists 
exclusively of up spins. 

The vectors can be expanded in the tensor product basis by the following pro¬ 
cedure: each arc is replaced by the spin singlet ||) (g) |J,) — |J,) (g) ||), and the unpaired 
spins are replaced by their tensor products. Finally, one symmetrizes in each block. 
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Fig. 7. Construction of a basis vector from an ordered Ising configuration 


Next, we briefly sketch how the properties claimed in Proposition 1 can be ver¬ 
ified. To do this we have to calculate the action of the Hamiltonian on the highest 
weight vector constructed in the previous paragraph. This is most easily accom¬ 
plished by deriving a graphical representation for the action of each term in the 
Hamiltonian as we did in the case of the pure spin 1 /2 chain. The Heisenberg inter¬ 
action for arbitrary spins of magnitude Sx and can be realized as an interaction 
between spin i ’s making up the spin Sx and Sa;+i, conjugated with the projections 
onto the symmetric vectors. The result is the following: 


—h. 


X,X-\-\ 




2 SxS- 


-Sx-S. 


xSx-{-l 


x + 1 


- 1 ) 



Here, the rectangles with label 2s represent the symmetrizing projections on the 
space of 2s spin | variables. The fundamental algebraic property that allows us to 
calculate the matrix elements of Hk graphically is the Jones-Wenzl relation (c.f., 
[10] and references therein): 



-k 


2s 

2s-k 1 



For any element of the basis introduced above one can now compute the action 
of the Hamiltonian and write it as a linear combination of the same basis vectors. 
From the grahical rules it is easy to observe that all off-diagonal matrix elements 
are non-positive. 

As before, it is straighforward to identify the basis for with a subset of the 
basis for The label of the rightmost box in any basis vector for the system 

k is raised by one but the number of arcs remains unchanged. 

The crucial property that allows us to compare the two Hamiltonians is the 
following. When Hk+i acts on a basis vector obtained from a corresponding Hk 
vector as we have just described, the only possible new terms that are generated 
are off-diagonal terms, which do not contain a bubble and, hence, are negative. 
The details of the calculation of these matrix elements and further applications will 
appear elsewhere [18]. 
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4 Extensions 

A highly desirable extension of our main result, of course, would be the proof of 
Conjecture 1 for all ferromagnetic isotropic Heisenberg models on an arbitrary graph! 
While we have been able to prove some partial results for the spin-1/2 model on 
an arbitrary tree and a few other graphs, we do not have an argument that works 
for arbitrary graphs [18]. But there are a few other directions in which one might 
extend the ordering of energy levels property. The aim of this section is to discus 
two such generalizations. In the first, the group SU{2) is replaced by the quantum 
group SUq{2), 0 < g < 1. This only works on the chain and, as far as we are aware, 
leads to information about physically interesting models in the case of the spin 1/2 
chain, namely the XXZ chain. The second generalization we consider is isotropic 
ferromagnetic models with higher order nearest neighbor interaction terms, such as 
{Sx ■ For this to be relevant, the spins have to be of magnitude > 1. 

4.1 The spin 1/2 S'Uq(2)-symmetric XXZ chain 

It is well-known that the translation invariant spin-1/2 XXZ chain with a particular 
choice of boundary helds is SUq{2) invariant [21]. This SUq{2) symmetry can be 
exploited in much the same way as the SU{2) symmetry of the isotropic model 
[12, 19]. Here we will show how it leads to a natural SUq{2) analogue of the FOEL 
property. 

The Hamiltonian of the S'[/q(2)-invariant ferromagnetic spin-1/2 chain of length 
L > 2 is given by 

L-l 

Hl = - + slsl+^) + {Slsl+i - 1/4)] (11) 

X = 1 

-A{A){Sl-S!). 

where Z\ > 1, and 

This model commutes with one of the two natural representation of SUq{2) on 
(C^)^, with q £ (0,1), such that A — {q + q~^)/2. Concretely, this means that Hl 
commutes with the three generators of this representation defined as follows: 

S® = ^ 111 (g) 

X = 1 

L 

x=l 
L 

S'-111® 

x=l 

where 



The SUq{2) commutation relations are 


• • • ® S® ® tx+i O ■ • • In 

• • • O tx-i ® S+ ® tx+i O ■ • • In 
■ ■ ■ 0 Sx ® ix+i ® • • • in 
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25 “^ — 25 “^ 

Note that one recovers the SU{2) definitions and commutation relations in the limit 
q —> 1. Hl also commutes with the Casimir opeator for SUq{2), given by 

c = s+s- + r = 

(g-i - g)2 


The eigenvalues of C are 

g-(2S+l) g2S+l 

(g-i - g)2 


S = 0,1/2,1,3/2,... 


and play the same role as S for the XXX model, e.g., they label the irreducible 
representations of SUq{2). The eigenspaces of C are invariant subspaces of and, 
as before, we denote the smallest eigenvalues of Hl restricted to these invariant 
subspaces by E{Hl, S). Note that the subspaces depend on q, but their dimensions 
are constant for 0 < g < 1. 


Theorem 4. 


E{Hl, S' + 1) < E{Hl, S), for all S < L/2 - 1. 

The proof of this theorem is identical to the one for the isotropic spin-1/2 chain 
up to substitution of the singlet vector ^ by the SUq{2) singlet fq = q |-l-) ® | —) — 

| —) ® |+), and changing the scalar value of the “bubble” to —{q + q~^). The details 
are given in [17]. 

4.2 Higher order interactions 

For spins of magnitude greater than 1/2 the Heisenberg interaction is not the only 
SU{2) invariant nearest neighbor interactions. It is easy to show that the most 
general SU{2) invariant interaction of two spins of magnitudes si and S 2 , i.e., any 
hermitian matrix commuting with the representation 0 is an arbitrary 

polynomial of degree < 2min{si,S2} in the Heisenberg interaction with real coeffi¬ 
cients: 

2 min{si ,S 2 } 

/112= ^ -52)"*. (12) 

m=0 

The definition of the FOEL property only uses S'[/(2)-invariance and therefore ap¬ 
plies directly to any Hamiltonian for a quantum spin system on a graph with at 
each edge an interactionof the form (12). We believe it is possible to determine the 
exact range of coupling constants such that FOEL holds for spin s chains with 
translation invariant interactions. So far, we have carried this out only for the spin-1 
chain. 

Theorem 5. FOEL holds for the spin-1 chains with Hamiltonian 

L-l 

X = 1 

with 0 < f3 < 1/3. Level crossings occur at fd = 1/3 and FOEL does not hold, in 
general, for (5 > 1/3. 
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The overall method of proof is the same as for the standard Heisenberg model. 
Theorem 2 applies directly, since only the SU{2) symmetry is used in its proof. The 
only difference is in the proof of Proposition 1. The same basis for the highest weight 
spaces is used but verifying the signs of the matrix elements is more involved. 


5 Applications 

In this section we discuss a number of results that are either consequences of the 
FOEL property, or other applications of the properties of the Heisenberg Hamilto¬ 
nian that allowed us to prove FOEL. 

5.1 Diagonalization at low energy 

The most direct applications of the FOEL property are its implications for the low- 
lying spectrum of the Hamiltonian. FOEL with strict inequality implies that the 
ground states are the multiplet of maximal spin which, of course, is not a new result. 
Since the maximal spin multiplet is unique, the first excited state must belong to less 
than maximal spin and therefore, by FOEL, to Smax — 1. In the case of models for 
which this second eigenvalue can be computed, such as translation invariant models 
on a lattice, this is again consistent with a well-known fact, namely that the lowest 
excitation are simple spinwaves. But in the case of arbitrary coupling constants and 
spin magnitudes in one dimension it proves that the first excited state is represented 
in the subspace of “one overturned spin” (with respect to the fully polarized ground 
state), i.e.. S'® = S^^x “ li which is a new result. 

More generally, the FOEL property can help with determining the spectrum of 
the Heisenberg model at low energies, whether by numerical or other means, in the 
following way. Suppose H is & Hamiltonian with the FOEL property. Diagonalize H 
in the subspaces for n = 0,1,... Ai, and select those eigenvalues that are 

less or equal than E{H, Smax — N). It is easy to see that the FOEL property implies 
that this way you have obtained all eigenvalues of the full H below < E{H, Smax—N). 

This is interesting because you only had to diagonalize the Hamiltonian in an 
invariant subspace that is explicitly known (by the representation theory of SU{2)) 
and of relatively low dimension: dim(Lfis 0{L"), while the full Hilbert 
space has dimension (2J -|- 1)® for L spin J variables. 

5.2 The ground states of fixed magnetization for the XXZ chain 

The spin 1/2 XXZ ferromagnetic chain with suitable boundary conditions, or de¬ 
fined on the appropriate infinite-chain Hilbert space has low-energy states that can 
be interpreted as well-defined magnetic domains in a background of opposite mag¬ 
netization [19, 11]. Using the techniques we have used for proving FOEL, we can 
rigorously determine the dispersion relation of a finite droplet of arbitrary size. 

The spin 1/2 XXZ chain can, in principle, be diagonalized using the Bethe Ansatz 
[13]. There are two complications that may prevent one from obtaining the desired 
information about its spectrum. The first is that a complete proof of completeness 
of the Bethe Ansatz eigenstates has been obtained and published only for the XXX 
chain {q — A = 1), although the corresponding result for the XXZ chain has been 
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announced quite some time ago [7]. The second problem is that the eigenvalues are 
the solutions of complicated sets of equations, such that proving statements as the 
one we discuss here, may be very hard. 

For brevity, let us consider the XXZ Hamiltonian for the inifinite chain defined 
on the Hilbert space generated by the orthonormal set of vectors representing n 
down spins in an infinite “sea” of up spins, and let us denote this space by Tin- 
Define 

E{n) = inf spec(i7|H„)- 

As before, the relation between q £ (0,1) and the anisotropy parameter A in the 
XXZ Hamiltonian (11) is given \yy A — {q + q~^)l‘l. 


Theorem 6. For n > I, we have 


(1 +( 72)(1 + 


Moreover, E{n) belongs to the continuous spectrum and is the bottom of a band of 


width 


4g" 


1-r 


(l + g")(l-g")' 


The states corresponding to this band can be interpreted as a droplet of size 
n with a definite momentum. The formula for the width indicates that the “mass” 
of a droplet diverges as n ^ oo. The proof of this result will appear in a separate 
paper [16]. If one looks back at the finite-volume eigenvalues E{Hk,k/2 — n) that 
converge to the bottom of the band in the infinite-volume limit {k oo, n fixed), the 
property (iv) amounts to property that the ground state energy of a droplet of fixed 
size n is strictly monotone decreasing in the volume. Moreover the finite-volume 
eigenvalues can be related to E(n), in the above limit, by using the generalization 
of the Perron-Frobenius result stated in Theorem 3. 


5.3 Aldous’ Conjecture for the Symmetric Simple Exclusion 
Process 

The Symmetric Simple Exclusion Process (SSEP) is a Markov process defined on 
particle configurations on a finite graph A. For our purposes it is convenient to define 
the process as a semigroup on Ha — where f^A is the space of configurations 

rj : A —f {0,1}. One thinks of ??(*) = 1 to indicate the presence of a particle at the 
vertex x. Let L be defined in Ha by the formula 

{Lf){'n) = Y^Y^Jxy{f{ri) - f{rt^^)) (13) 

r) 

where denotes the configuration obtained form q by interchanging the values of 
q{x) and q{y). The parameters Jxy are positive numbers representing the jump rate 
at the edge x ^ y € A. 

Clearly, the number of particles is a conserved quantity of the process. Con¬ 
cretely, this means that Ha decomposes into a direct sum of invariant subspaces 
H^a \ n = 0,..., |A|, where 77^"^ consists of all functions supported on configura¬ 
tions q that have extactly n particles, i.e., ’7(®) = n. 0 is a simple eigenvalue of 
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each of the restrictions L\ („), and L is non-negative dehnite. For each n, the corre- 
spending invariant measure is the uniform distribution on n-particle configurations. 
Let A(n) denote the smallest positive eigenvalue of L\ in). Since the dynamics 

of the SSEP is given by the semigroup {e“*^}t>o, A(n), determines the speed of 
relaxation to the invariant measure. 

The following conjecture is known as Aldous’ Conjecture but his website [1] 
Aldous states that it arose in a conversation with Diaconis. So, maybe it should be 
called the Aldous-Diaconis Conjecture. 

Conjecture 2. 

A(n) = A(l), for all 1 < n < |A| — 1. 

Apart from being a striking property, namely that the relaxation rate should be 
independent of the number of particles, it could also be very useful. The SSEP for 
one particle is just a random walk on the graph A, and many powerful techniques 
are available to study the relaxation rate of random walks. In physical terms one 
would say that the conjectured property reduces the many-body problem of finding 
the relaxation rate for n particles to a single one-particle problem. 

Proposition 2. If the ferromagnetic spin-1/2 Heisenberg model with coupling con¬ 
stants Jxy on a graph A satisfies FOEL, then Conjecture 2 holds for the SSEP on 
A with jump rates Jxy 12.. 

Proof. The proof is based on the unitary equivalence of L and the ferromagnetic 
spin 1/2 Heisenberg Hamiltonian H given by 

x^yeA ■' 

The unitary transformation U : LP'{Qa) —> Ha ~ that relates L and H is 

explicitly given by 

9 f r^Uf = ip = Y fiv) \v), where Sl \rj) = (rjx - 1/2) \rj) 

V 

To see this note that 

1/4-5, • Sy = {l-txy)/2, 

where txy interchanges the states at x and y in any tensor product vector. Then 

"v-'jEE f{v)Jxy{l - txy) \V) 
xc-^y r) 

- firT^))\v) 

xc-^y 7) 

v 

Therefore, HUf = ULf, for all / G L'^{nA. 

Under this unitary transformation, the particle number becomes the third com¬ 
ponent of the total spin: 
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S'tot — ~ 1^1/2 + n. 

The unique invariant measure of SSEP for n particles is the uniform measure on 
{ri £ Ha \ The corresponding state for the spin model belongs to 

the unique multiplet of maximal total spin, i.e., is a ground state. A(n) is the next 
eigenvalue of H is the same value of total S'®. Since the first excited state of H, 
by FOEL, is a multiplet of total spin Smax — 1, this eigenvalue has an eigenvalue 
with any value of S® in the range —Smax + 1,..., Smax — 1. We have Smax = |T|/2. 
Therefore, this corresponds to the range 1 < n < |/1| — 1. Hence, A(n) is independent 
of n in this range. 

In combination with Theorem 1, this proposition has the following corollary. 
Corollary 1. Conjecture 2 holds for ehains. 

Our partial result for trees (not discussed here) also implies Conjecture 2 for 
arbitrary finite trees as well as some graphs derived from trees. These cases of 
the Aldous-Diaconis conjecture were previously know [2, 8], as well as some other 
examples where one can compute A(n) exactly [3, 5, 4]. Needless to say, a full proof of 
FOEL, the Aldous-Diaconis Conjecture, or even a proof for additionial special cases, 
would be of great interest. An interesting direction for generalization considered by 
Aldous is to also establish the analogous formula for the spectral gap for a card- 
shuffling model with full SU{n) symmetry, which restricts to the SSEP when one 
considers cards of only two colors. 
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